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1. Introduction 

'Rare decays' of B mesons, such as, B ^ K*y, B ^ K*}X^ji.^, B ^ py, ... where ^ ^ 57 are 
flavour changing neutral current or FCNC processes and are thus not aflowed at tree level by the 
GIM mechanism. However this makes them sensitive to higher scales, and may affect various CKM 
matrix elements, such as Vts or Vtd- These exclusive events can be investigated at the LHC by the 
LHCb experiment. A theoretical framework is provided by QCD factorisation, eg, (which 
is a heavy quark expansion in I /mi,), perturbative QCD [^, soft-coUinear effective theory [Q] or 
light-cone sum rules [||]. These give a decay amplitude related to vector distribution amplitudes or 
vector DAs. These are usually defined in the MS scheme at some scale /i. In this article we compute 
using lattice QCD the lowest moment of the K* DA. Analogous computations have recently been 
performed for the spin particles n and K, [§, ^. 

As we have vector particles, with a polarisation vector, we have two distinct DAs: 
and 0^(<^). These are functions of i§ G [— 1,+1], where x = \{l + t,) and I —x = ^(l — ^) are 
the fractions of meson momentum carried by the quark and anti-quark respectively (in the infinite 
momentum frame). An expansion in terms of Gegenbauer polynomials 



with 



allows (possible) reconstruction of the full DA. In particular as ajl'^ when jj. —K>a, v/e might 
hope that knowledge of the lowest lowest few aj!'^ coefficients suffices. Indeed the lattice compu- 
tation is only capable of giving low moments of DAs, defined by 

where a\"^ = 5/3(i^)ll'^, a^'^ = 7/12(5(<^^)ll'-'- — 1), . . . .As Gegenbauer polynomials are orthog- 
onal polynomials with weight l — ^^ and as Cq^^ = 1 then the normalisation is such that (1)"'^ = 1. 
Finally we note that G-parity restricts the functional form of (j)p'^ to an even function of and so 
non-zero moments are {^)k*, {^^)k*> {^^)p' •••• 

2. Minkowski matrix elements 



Longitudinal matrix elements are given by 

.11 



with 



2 
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where q = d or s, means symmetrised and traceless in these indices, ^= 
polarisation index. Correspondingly transverse matrix elements are given by 



and A is the 



(where a^'^^^" = \ [y<-«')M^y(-^)v]) ^j^h operators 



This all looks rather complicated, but for no derivatives {n = 0) the equations reduce to the familar 
ones for the Fy^ decay constants, namely 



X ■ 



V 



and 



Thus we see that these equations have been normalised with Fy to ensure, as required, that 

(1)11.^ = 1. 



3. The Lattice 

On the lattice we need a careful choice of lattice operators to avoid mixing with same di- 
mension operators, and worse mixing with lower dimensional operators when I /a subtractions are 
required. We shall consider only n = 1 operators here, the list ^ used is 



n 


Operator 


Representation 


1 




Tr,C = + 


1 







for the longitudinal operators, where 



n 


Operator 


Rep. 


1 


fff-f" = 6 




^Aji, i + j 


^2 ' 
^1 ' 


c = 


+ 


1 


^^■^ = t 






c = 


+ 



for the transverse operators, where ffv^lo^l^■■■^„ = ^7v7mo ^mi^M2 " ' " ^^^„ " / Mo)- The operators 
belonging to different (hypercubical) representations have been labelled by 'a' and 'b', and should 
give the same results, at least in the continuum limit. (Further results, including n = 2 operators 
will appear in [^.) 

Correlation functions are then defined, where 



cMt;p) = m;p)aiO;p)'^), 
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with Q. = V orT, where to improve the signal these operators have been 'Jacobi' smeared. Then in- 
serting complete sets of states in the standard way gives correlation functions involving {0\^\V,p,X) 
and {0\ff\V,p,X). The unwanted {0\Q.\V,p,X) may be cancelled by forming ratios. For example 
we find for some of the (bare) operators 

• Longitudinal 



3i5lnPm{^)a 

-Ev (^^1^) coth£v(^^, -0 

and similar expressions for the other operators as (in the above V) can also be replaced by T 
giving further ratios. Thus many cross checks are possible. Note that the t fit function is known and 
may be either tanh, coth or 1. Also half the n = I operators can be measured at zero momentum; 
the others cannot. However for those operators a non-zero ratio requires only a single unit of 
momentum in one direction. We choose the lowest possible momentum, \p\ = 2%/Ns and average 
over the three spatial directions. 

We use unquenched «/ = 2, 0{a) improved clover fermions in our simulations, the lattices 
employed being: 



i8 


^sea 


nIxNt 


Trajs 


nips /my 




a[fm] 




mps[MeV] 


5.29 


0.1350 


16^ X 32 


5700 


0.76 


6.7 


0.075 


1.20 


1100 


5.29 


0.1355 


24^ X 48 


2100 


0.70 


7.8 


0.075 


1.81 


860 


5.29 


0.1359 


24^ X 48 


4900 


0.62 


5.7 


0.075 


1.81 


630 



7LiC^\\.ay{t;p) 

jUiCv.v.it'P) 
C^lbv_{t;p) 

j'LiCv,v,{t;p) 

• Transverse 

^HiC ^±,byit; p) 



together with various Ky^i for the valence quarks. Note that L5 = aNs and mj^+ /mp+ ~ 0.18. The 
scale is set from the ro force scale, using a value of of ro = 0.467 fm = 1/422.5 MeV. a is deter- 
mined from extrapolating (ro/a) to the chiral limit (presendy giving {ro/a)c{P = 5.29) = 6.20(3)). 
No operator improvement has been attempted, although experience from quenched unpolarised op- 
erators has indicated that these effects are probably small, [10|. 

A non-perturbative renormalisation - RI' — MOM method has been used to determine the 
renormalisation constants. {Z'^''^°" is computed numerically and from this Z'^' is determined. This 
is then converted to Z'^^{ix = 2GeV), which is the scheme and scale that all our results are presented 
here.) For more details see the forthcoming paper [|lT]]. 

A (typical) result for the ratio is shown in Fig. [T], where we observe a clear tanh function. 
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Figure 1: The ratio -;|;|l;C^°v,(f;0)/ii:,Cy,v,(f;6) = {^)ltanhmv{jNT - 1) versus f for j3 = 5.29, 
T^sea = 0.1350, Kyai = (0.1355,0.1343), p — 0. The results are denoted by filled circles. Also shown is a 
one parameter fit (the my mass having been determined previously). The fit range is denoted by vertical 
dashed lines. The operator has been renormalised to the MS scheme at /i = 2 Ge V, so that the value obtained 
corresponds directly to a point in Fig. ||(the sixth point from the left). 

4. Results 

As noted previously, odd moments vanish for degenerate mass (valence) quarks and thus we 
have (m^2 O^gi) 

oc mlp^-ml^, 

where nips is a pseudoscalar meson with degenerate mass quarks and rriKps is a pseudoscalar meson 
with possibly non-degenerate mass quarks. (For the even moments, not considered here, there is 
no such restriction and are just symmetric in the quark masses.) For (<^)|-»^ we first, for fixed m^ga, 
plot (<§)"'^ against (valence pseudoscalar masses) — m^^ and interpolate to the physical point 
mj^ — m|, This is then taken as a function of ntsea °= m^s ^"^^ extrapolated to the chiral Umit to 
give finally 

In Fig. ^ we show versus m^p^ — m^^ together with a one-parameter fit passing through 
the origin. Also shown (red star) is the value when mj^^^ — m^^. = mj. — m\. Fig. ^ shows the 
corresponding results for 

As discussed previously we must now extrapolate nisea °^ ni^ps to the chiral limit (the difference 
between this and m| is negligible). In Fig. ^ we show this extrapolation for giving an estimate 
for In Fig. || we show the equivalent picture for (i^)j- leading to a value for (<^)^*. 

This is repeated for other channels and we thus finally arrive at the (preliminary) results 

« 0.033(2)(4) ajj,, « 0.055(3)(7) 

i^)-^, w 0.030(2)(8) af^, ^ 0.050(3)(13) 
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Figure 2: (^)I versus mj.^^^ - m^j in the MS scheme at /i = 2GeV for j3 = 5.9, K^ea = 0.1350, p = for 
various valence quark combinations. A linear fit vanishing when the two valence quark masses are the same 

2 2 2 2 

Kps ~ '^ps ^ '^K ~ 



is also shown. The red star shows the value when "-^ ^ '^^ 
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Figure 3: (^)J^ versus m^^^ — m^,. Same notation as for Fig. ^. Note that here we work at finite momentum 
\p\ = 2n/Ns. 
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Figure 4: (^)J| versus m^, for the three sea quark masses K^ea = 0.1350, 0.1355 and 0.1359 (black circles), 
together with a linear extrapolation to the chiral limit (red circle). 



0.04 r 




Figure 5: {t,)-^ versus rni^. Same notation as in Fig. 



6 



Distribution Amplitudes of Vector Mesons 



R. Horsley 



(in the MS'-scheme at a scale of = 2GeV) where the first error comes from the spread of channels 
presently analysed and the second error is an estimation of possible chiral extrapolation error (the 
fit being repeated dropping one data point). Also any discretisation errors have been ignored. 

These are to be compared with the results from sum rule estimates of 0.02(2), aj^/^, 

0.03(3) [ jl2| ] at the same scale, and the limit function ^"•'^y"'P(^) giving a\'^, = 0. Potentially lattice 
results are more reliable than sum rule estimates and may help in a reconstruction of the vector 
distribution amplitude. 

Our conclusion is that a lattice determination of (moments of) vector DAs is possible. We plan 
to extend these results to lighter pseudoscalar masses, j3 = 5.40 (a finer lattice) and to (i^^)"'^ for 
both the K* and p. Further results (including the zero moment decay constants) will appear in [Q]. 
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